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Abstract. A non-perturbative approach to Generalized Parton Distribu- 
tions, and to Deeply Virtual Compton Scattering in particular, based on 
off mass shell extension of dual amplitudes with Mandelstam analyticity 
(DAMA) is developed with the spin and helicity structure as well as the 
threshold behavior accounted for. The model is tested against the data on 
deep inelastic electron- proton scattering from JLab. 
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1. Introduction 

Parton distributions measure the probability that a quark or gluon carry a 
fraction x of the hadron momentum. Relevant structure functions (SF) are 
related by unitarity to the imaginary part of the forward Compton scatter- 
ing amplitude. Generalized parton distributions (GPD) [1, 2, 3] represent 
the interference of different wave functions: one, where a parton carries mo- 
mentum fraction x + ^, and the other one, with momentum fraction x — ^, 
correspondingly; here ^ is skewedness and can be determined by the exter- 
nal momenta in a deeply virtual Compton scattering (DVCS) experiment. 
Basically, the DVCS amplitude can be viewed as a binary hadronic scat- 
tering amplitude continued off mass shell. DVCS is related to the GPD in 
the same way as elastic Compton scattering is related to the ordinary SFs. 

Apart from the longitudinal momentum fraction of gluons and quarks, 
GPDs contain also information about their transverse location, given by 
a Fourier transformation over the t-dependent GPD. Real-space images 
of the target can thus be obtained in a completely new way [4]. Spacial 
resolution is determined by the virtuality of the incoming photon. Quantum 
photographs of the nucleons and nuclei with resolutions on the scale of a 
fraction of femtometer are thus feasible. 

One of the first experimental observations of DVCS was based on the 

recent analysis of the JLab data from the CLAS collaboration. New mea- 
surements at higher energies are currently being analyzed, and dedicated 
experiments are planned [5]. 

On the theoretical side, much progress has been achieved [1, 2, 3] in 
treating GPD in the framework of quantum chromodynamics (QCD) and 
the light-cone technique. On the other hand, non-nonpertubative effects 
(resonance production, the background, low-Q^ effects) dominating the 
kincmatical region of present measurements and the underlying dynam- 
ics still leave much ambiguity in the above-mentioned field-theoretical ap- 
proach. Therefore, as an alternative or complementary approach we have 
suggested [6, 7, 8] to use dual amplitudes with Mandelstam analyticity 
(DAMA) as a model for GPD in general and DVCS in particular. We remind 
that DAMA realizes duality between direct-channel resonances and high- 
energy Reggc behavior ( "Veneziano-duality" ) . By introducing Q^-dependence 
in DAMA, wc have extended the model off mass shell and have shown [6, 7] 
how parton-hadron (or "Bloom-Gilman" ) duality is realized in this way. 
With the above specification, DAMA can serve as and explicit model valid 
for all values of the Mandelstam variables s, t and u as well as any Q^, thus 
realizing the ideas of DVCS and related GPDs. The historical and logical 
connections between different kinametical regions and relevant models are 
depicted in Fig. 1. 
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Figure 1. Generalized parton distribution amplitudes. A road map. 



In this work we concentrate on the very delicate and disputable problem 
of the off mass shell continuation (introduction of the variable Q^) in the 
dual model, starting with inclusive electron-nucleon scattering both at high 
energies, typical of HERA, and low energies, with the JLab data in mind 
(see ref. [8] for more details). 

2. Simplified model 

The central object of the present study is the nucleon SF, uniquely related 
to the photoproduction cross section by 

47ra(l + 

where the total cross section, a] ^, is the imaginary part of the forward 
Compton scattering amplitude, A{s,Q'^), a] ^(s) =1m A(s,(5^); m is the 
nucleon mass, a is the fine structure constant. The center of mass energy of 
the 7*p system, the negative squared photon virtuality and the Bjorken 
variable x are related hy s = Q'^{1 — x)/x + m?. 

We adopt the two-component picture of strong interactions [9] , accord- 
ing to which direct-channel resonances are dual to cross-channel Regge 
exchanges and the smooth background in the s— channel is dual to the 
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Pomeron exchange in the i— channel. This nice idea, formulated [9] more 
than three decades ago, was first realized explicitly in the framework of 
DAMA with nonlinear trajectories (see [11, 6] and earlier references therein). 

As explained in Refs. [6] and [11], the background in a dual model 
corresponds to a pole term with an exotic trajectory that does not produce 
any resonance. 

In the dual-Rcgge approach [6, 7, 8] Compton scattering can be viewed 
as an off mass shell continuation of a hadronic reaction, dominated in the 
resonance region by direct-channel non-strange {N and A) baryon trajec- 
tories. The scattering amplitude follows from the pole decomposition of a 
dual amplitude [6] 



Ais,Q^) 



norm ^ A E. ' 

mm *■ V / 



t=0 i=N*,N*,A,E 



where i runs over all the trajectories allowed by quantum number exchange, 
norm and Ai^s are constants, fiiQ'^Ys are the form factors. These form 
factors generalize the concept of inelastic (transition) form factors to the 
case of continuous spin, represented by the direct-channel trajectories. The 
j^min j.gfgj.g lY^Q gpj]^ of ^]jg £j.g^ resonance on the corresponding trajectory 
i (it is convenient to shift the trajectories by 1/2, therefore we use aj = 
^phys which due to the semi-integer values of the baryon spin leaves 

n in Eq. (2) integer). The sum over n runs with step 2 (in order to conserve 
parity). 



It follows from Eq. (2) that 



[/;(Q^)]^^"-"'""+^^X" 
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Xm ^(s, (5^) = norm ^ E 

i=N*mA.E n=n^i^ ~ «i(s))^ + (Xm ai{s)) 

(3) 

The first three terms in (3) are the non-singlet, or Reggeon contributions 
with the A^* and A trajectories in the s-channel, dual to the exchange of an 
effective bosonic trajectory (essentially, /) in the t-channel, and the fourth 
term is the contribution from the smooth background, modeled by a non- 
resonance pole term with an exotic trajectory a£;(s), dual to the Pomeron 
(see Ref. [6]). As argued in Ref. [6], only a limited number, M, of resonances 
appear on the trajectories, for which reason we tentatively set J\f = 3 - one 
resonance on each trajectories (A^f , A|, A), i.e. n™"^' = n""". Our analyses 
[8] shows that = 3 is a reasonable approximation. The limited (small) 
number of resonances contributing to the cross section results not only 
from the termination of resonances on a trajectory but even more due to 
the strong suppression coming from the numerator (increasing powers of 
the form factors). 
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TABLE 1. Values of the fitted parameters. In the first column 
we show the result of the fit when the parameters of the bary- 
onic trajectories are fixed. The second column contains the result 
of the fit when the parameters of the trajectories are varied. ^ de- 
notes the parameters of the physical baryon trajectories from ref. [8]. 

* The coefficient norm is chosen in such a way as to keep 
^iv* = 1 in order to sec the interplay between different resonances. 

* Using intercepts and thresholds as a free parameters does not im- 
prove the fit, but they may get values far from original. 







-0.8377 (fixed)'f 


-0.8070 


-0.8377 (fixed)* 




Ql 
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0.85 (fixed) t 
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A 




0.1969 (fixed)t 
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1.5722 
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0.5484 
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0.1126 
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2.691 




Q'o' 
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norm 


0.021 


0.0207 


0.0659 




Xd.o.f. 


28.29 


11.60 


18.1 



We use Regge trajectories with a threshold singularity and nonvanishing 
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imaginary part in the form: 

a{s) = ao + ais + a2{^/s^ - s/s^^^), (4) 

where sq is the hghtest threshold, sq = ("J'Tr+'Tip)^ = 1.14 GeV^ in our case, 
and the hnear term approximates the contribution from heavy thresholds 
[6, 7, 8]. 

For the exotic trajectory we also keep only one term in the sum ^. n^*" 
is the first integer larger then Max{TZe oe) - to make sure that there are 
no resonances on the exotic trajectory. The exotic trajectory is taken in the 
form 

aE{s) = aE{0) + aiE{^/sE - V^e- s) , (5) 

where Q!e(0), aiE and the effective exotic threshold se are free parameters. 

To start with, wc use the simplest, dipole model for the form factors, 
disregarding the spin structure of the amplitude and the difference between 
electric and magnetic form factors: 

MQ')={l + QVQl,y' , (6) 

where Qq ^ are scaling parameters, determining the relative growth of the 
resonance peaks and background. 

We test our model against the experimental data from SLAG [12] and 
JLab [15]^. This set of the experimental data is not homogeneous, i.e. 
points at low s (high x) are given with very small experimental errors, 
thus "weighting" the fitting procedure not uniformly. Consequently, a pre- 
selection of the date involved in the fitting procedure was made, although 
we present all the experimental points in the Figure. The results of our fits 
and the values of the fitted parameters are presented in Figs. 2 and Table 
1. More details can be found in ref. [8]. 

For the first fit - first column of Table 1 and dashed-dottcd lines in 
Fig. 2 - we fix the parameters of the A and A^* trajectories to the physical 
ones, reproducing the correct masses and widths of the resonances, leaving 
the four scaling constants Qf, four factors Ai and the parameters of the 
exotic trajectories to be fitted to the data. Then, to improve the model we 
effectively account for the large number of overlapping resonances (about 
20) present in the energy range under investigation. For this reason we con- 
sider the dominant resonances {N^, N2 and A) as "effective" contributions 
to the SF. In other words, we require that they mimic the contribution of 

^In Ref. [11] the whole DAMA integral was calculated numerical and it has been shown 
that in the resonance region the direct-channel exotic trajectory gives a non-negligible 
contribution, amounting to about 10-12%. 

^We are grateful to M.I. Niculescu for making her data compilation available to us. 
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the dominant resonances plus the large number of subleading contributions, 
which, together, fully describe the real physical system. In the light of these 
considerations, we have refitted the data, allowing the baryon trajectories 
parameters to vary. The resulting parameters of such a fit are reported in 
Table 1 (second column) . It is worth noting that although the range of vari- 
ation was not restricted, the new parameters of the trajectories stay close 
to their physical values, showing stability of the fit and thus reinforcing our 
previous considerations. From the relevant plots, shown in Fig. 2 with full 
lines, one can see that the improvement is significant, although agreement 
is still far from being perfect {Xdof ~ H-S)- 

The smooth dashed lines in Fig. (2) correspond to a "scaling curves" , i.e. 
a phenomenological parameterizations of the SF exhibiting Bjorken scaling 
and fitting the data [10]. 

3. Spin and helicity; threshold behaviour 

In the preceding section the nuclcon form factor was treated in a phe- 
nomenological way, fitted to the data, neglecting the longitudinal cross 
section ctl, and the ^ limit of the SFs. On the other hand, it is 
known [14, 13] that with account for both ai and ctt the form factor can 
be presented as sum o three terms: G+(Q^). Go(Q^) and Go(Q^), corre- 
sponding to 7*A^ R helicity transition amplitudes in the rest frame of 
the resonance R: 

Ga^ = ^ < R, Xr = \n- A^I J(0)|iV, Xn >; (7) 

here Xr, \n and A-y are the resonance, nucleon and photon helicities, J(0) is 

the current operator; A^ assumes the values —1, and -|-1. Correspondingly, 
one replaces the Q^— dependent expression in the numerator of (3) by 

fiQY ^ |G+(g2)|2 + 2|Go(Q2)|2 + |G_(Q2)|2. (8) 

The explicit form of these form factors is known only near their thresh- 
olds |g| ^ 0, while their large-Q^ behavior may be constrained by the quark 
counting rules. 

According to [16], one has near the threshold 

g±(q2) ~ Go(Q2)~|i,r-v2 (9) 

for the so-called normal (1/2+ — > 3/2~, 5/2+, 7/2~, ...) transitions and 



(10) 
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for the anomalous (1/2^ 1/2 ,3/2"*", 5/2 ,...) transitions, where 

,^ v/(M2 - - Q2)2 + 4M2Q2 ,^ - 

k1 = 2M ' l^"1'3=° = ' 

M is a resonance mass ^. 

FoUowing the quark counting rules, in refs. [13] (for a recent treatment 
see [14]), the large-Q^ behavior of G's was assumed to be 

G+(g2)^g-3, Go(Q2)~Q-^ G'_(Q2)^q-5. (13) 

Let us note that while this is reasonable (modulo logarithmic factors) for 
elastic form factors, it may not be true any more for inelastic (transition) 
form factors. For example, dual models (see Eq. (1) and ref. [6]) predict 
powers of the form factors to increase with increasing excitation (resonance 
spin). This discrepancy can be resolved only experimentally, although a 
model-independent analysis of the Q^.^jepgndence for various nuclear exci- 
tations is biased by the (unknown) background. 

In ref. [14] the following expressions for the G's, combining the above 
threshold- (9), (10) with the asymptotic behavior (13), was suggested: 



' °' W + Ql \e[\Q\Q=oQ' + Q'i W + Ql 



for the normal transitions and 



i^^i^-'^^^^^Km^q^J MWJ '''' 



' °' W + Ql) [WQ' + Q'oV W + Ql) ^ ^ 

^In our approach M is defined by trajectory: TZea{M^) = J, where J is a spin of the 
resonance; and therefore for the varying baryon trajectory M might differ from the mass 
of physical resonance. 
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for the anomalous ones, where m+ = 3, nio = 4, m_ = 5 count the quarks, 
C and a are free parameters. The form factors at = are related to the 
known (measurable) helicity photoproduction amplitudes ^1/2 ^3/2 by 

l^-C)!"^/^!^'/../.!. (18) 

The values of the helicity amplitudes are quoted by experimentalists [17] 
(those relevant to the present discussion are compiled also in [14]). 

We have fitted the above model to the JLab data [15] again by keep- 
ing the contribution from three prominent resonances, namely A(1232), 
A''*(1520) and A''*(1680). For the sake of simplicity we neglected the cross 
term containing Go (and the coefficient C and parameter a), since it is 
small relative to the other two terms in Eq. (8) [14]. 

We use the same background term as in the previous section, but with 
the normalization coefficient Ge- 



background = . , , !e = Ge^ 



To be specific, we write explicitly the three resonance terms, to be fitted 
to the data, (cf. with Eq. (3) of the previous section)^: 

Tm A{s,Q^) = norm- 
f^^^ + fN-{N- + fN+lN^ ^ background, 



+ (l-E^_)2 + /2^_ (2-i2^+)2 + j2^^ 

(19) 

where e.g. /a is calculated according to Eqs. (8,14,16): 

(20) 

here R and I denote the real and the imaginary parts of the relevant tra- 
jectory, specified by the subscript. Similar expressions can be easily cast for 
/jv+ and for as well. 

The form factors at = can be simply calculated from Eq. (18) by 
inserting the known [17] (see also Table 1 in [14]) values of the relevant 
photoproduction amplitudes: 

^Remember that our trajectories are shifted by 1/2 from the physical trajectories 
ai = af' - 1/2. 
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^A(i232)(l/2,3/2) = (-0.141 -0.258 GeV-^l'')- 

^iV(i520))(l/2,3/2) = (-0.022 GeF-V2, 0.167 GeV'^l''); 

^Ar(i680)(l/2,3/2) = (-0.017 GeF-V2, 0.127 GeF-V2). 

Note that in this way the relative normahzation of the resonance terms 
is fixed (no Aj's appearing in previous section), leaving only two adjustable 
parameters, Q\ and Q'q (provided C is set zero, then a also disappears), 
which means that this version of the model (with the cross term, containing 
Go, neglected!) is very restrictive. 

The resulting fits to the SLAC and JLab data are presented in Fig. 2 
(dashed lines) and Table 1 (third column) . The fit is not so good (x^ o f ~ 
18.1), what probably tells us not to neglect the cross term. We hope to 
account for this and to improve the fit in the forthcoming works. 

4. Future Prospects 

In this work we have presented new results on the extension of a dual 
model (Sec. 2) to include the spin and helicity structure of the amplitudes 
(SFs) as well as its threshold behavior as — > 0. Let us remind that the 
lowest threshold (in the direct channel) of a hadronic reaction (if Compton 
scattering is considered in analogy with ttA^ scattering) is sq = (m7r+mjv)^. 
Hence its imaginary part and the relevant structure function starts at this 
value and vanishes below sq- The "gap" between the lowest electromagnetic 
TTT-p and the above-mentioned hadronic thresholds is filled by the arguments 
and formulas of Sec. 3. The — > is also important as the transition 
point between elastic and inelastic dynamics. 

The important step performed in our work after Ref. [14] is the "Reg- 
gization" of the Breit-Wigner pole terms (19), i.e. single resonance terms 
in Ref. [14] are replaced by those including relevant baryon trajectories^. 
The form of these trajectories, constrained by analyticity, unitarity and by 
the experimental data is crucial for the dynamic. The use of baryon trajec- 
tories instead of individual resonances not only makes the model economic 
(several resonances are replaced by one trajectory) but also helps in clas- 
sifying the resonances, by including the "right" ones and eliminating those 
nonexistent. The construction of feasible models of baryon trajectories, fit- 
ting the data on resonances and yet satisfying the theoretical bounds, is in 
progress. 

Note also that the dual model constrains the form factors: as seen from 
Eq. (3), the powers of the form factors rise with increasing spin of the 
excited state. This property of the model can and should be tested experi- 
mentally. 

^Another important step is adding the background into consideration. 
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The dependence, introduced in this model via the form factors can 
be studied and constrained also by means of the QCD evolution equations. 

To summarize, the model presented in this paper, can be used as a 
laboratory for testing various ideas of the analytical S-matrix and quantum 
chromodynamics. Its virtue is a simple explicit form, to be elaborated and 
constrained further both by theory and experimental tests. 

Acknowledgment. L.J., A.L. and V.M. acknowledge the support by IN- 
TAS, grants 00-00366 and 97-31726. 
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